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ABSTRACT 

A lineeu:  spline  method  for  the  solution  of  the  Abel  integral  equation 
X , 

/ , y(s)ds  « f (x)  , X ^ 0 

0 /x^-s^ 

is  analyzed.  The  approximate  solution  along  with  its  derivative  converges  to 
the  corresponding  exact  solutions  at  each  point  in  the  interval  of  integration, 
the  orders  of  convergence  being  two  and  one,  respectively.  An  asymptotic  for- 
mula for  the  discretization  error  is  obtained.  The  method  is  illustrated  by 
a numerical  example. 


AMS  (MOS)  Subject  Classifications;  45E10,  45L10,  65R05. 

Key  Words:  Abel  Integral  equation.  Global  approximation.  Linear  spline  method. 
Asymptotic  error  formula. 
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Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024 . 


Abel-type  integral  equations  (see  Abstract)  occur  regularly  in  applications.  i 

Typical  examples  are  the  determination  of  the  emission  coefficients  in  radiation  I 

technology,  the  determination  of  gravitational  anomalies  for  an  axially  symmetric 

distribution  of  masses,  and  the  analysis  of  the  fringe  shifts  in  interferograms . 

Applications  of  Abel  enuations  are  usually  directly  or  indirectly  related  to 

systems  with  axially  or  radially  symmetric  geometry.  i 

The  equation  considered  in  this  paper  has  an  explicit  inversion  formula.  i 

Numerical  methods  based  on  this  formula  have  been  investigated,  but  they  all  have  [ 

I 

to  deal  with  the  presence  of  a derivative  in  the  inversion  formula.  Experience  . 

shows  that  direct  methods  are  almost  as  effective  as  using  the  inversion  formula, 
and  they  can  be  generalized  to  solve  equations  for  which  an  inversion  formula 
is  not  known.  In  particular,  the  results  developed  in  this  paper  give  some  idea  | 

of  the  necessary  tools  and  possible  results  for  linear  spline  or  higher  degree  i 

spline  approximate  solution  for  a more  general  Abel  integral  equation  of  the 


form: 

X 

J ?LiJ y(s)ds  = f(x)  , 0<a<l,  0<_B^  1-a,  0 ^ x £ 1 

0 (x-s)™(x+s)^ 

The  trapezoidal  product  integration  method  considered  here  has  been  analyzed 
by  Atkinson  and  Benson.  The  contribution  of  the  present  paper  is  that  we  can  ob- 
tain similar  results  under  slightly  weaker  assumptions  by  a much  simpler  method. 

The  asymptotic  error  formula  obtained  for  the  numerical  solution  is  simple 
and  quite  remarkable  because,  if  we  are  solving  differential  or  integral  eaua- 
tions  numerically,  we  usually  have  to  solve  another  differential  or  integral 
equation  to  get  an  error  estimate,  whereas  for  Abel  equations  of  the  type  con- 


•plin*  function  P(x),  with  knots  at  tha  points  as  an  aporoximatlon  to  y{x) , i.o., 

for  I • 0,  I,  2,»*- 


(J.l> 


P(x)  - - - x)V^  t (X  - Xi)Yi^jl.  x^ixlx^^j  . 


Tha  function  P(x)  Is  continuous  at  the  knots. 

The  approximate  solution  of  the  integral  equation  Is  obtained  by  requiring  that  (1.1) 
be  satisfied  at  the  knots  x^,  i.e.,  the  exact  solution  y(x)  is  replaced  by  the  approximate 

solution  P(x)  derived  from  the  value  P(x^)  • computed  from 


(2.2) 


X , 

f — -=rP(s)ds  - f(x  ),  k “ 1,  2,“’ 

0/22 

/x^-s 


This  can  be  rewritten  in  the  form: 

k 


(2.3) 


“k.i  ^ “ *‘V' 

i«0 


k - 1,  2,- 


where 


(2.4) 


1 (X  -s) 

“k.o  “ L n^t  h ' 


■‘0 


{ “k.i  ■ I h 

1-1 


‘-■Vi’  . , 1 ‘*in-®’ 

ds  V I — • 7 


'i 


'2  2 


ds,  i “ 1,  • • ‘/k  - 1 , 


I 1 

“ L /T~2 


(s-Xk.i) 


ds 


\-l 


Equation  (2.1)  la  a nonsinqular  triangular  system  for  Y^,  since 


(2.5) 


4 1 4 1 

3 T^i-k.kii 


for  k “ 1,  2,***.  The  starting  value  of  this  system  needs  to  be  determined  by  other  means, 
for  instance,  it  might  be  obtained  from 


(2.6) 


Yq  - V(0)  - - f(0)  . 


-2- 


which  «xlitta,  wh«*nev»i  (1.1)  has  a continuous  solution. 

An  astimata  of  y'(x)  is  qivan  hy  tha  derivative  of  (2.1).  If  wa  denote  this  (constant) 
estimate  of  y'(x)  in  x^  ^ x < x^^j  bv  Y|,  this  qives 


(2.7) 


Y*  • — lY  - Y ) 
i h ' itl  ''i'* 


*i  - * ' ’‘i+1 


I . Oonvatqe nee  of  the  Method . 

The  proofs  of  Theorem  1.1  and  Theorem  4.1  require  the  followinq  lemma. 

**  there  exist  a ixinstant  C ' 0 ami  an  inteqer  N 2.  1 • «H  indenendent  of  )t, 

such  that 


' C. 


with 


t lien 


^ ‘ 1^1- 


‘ J'Vi.il  ' 

i«0 


1^1 


I • 0,  !,•••, N , 

k - N,  Nrl,**", 


k • N,  N* ■ 

i « 0,  1,2,- 


I'nxif : Assume  that 


I X j I ^ t'  for  i -O,  l,'*»,k  (k2.N)-  T)ien 
k k 

1^-'  F l'\*.  J - I J'  ■ ^ • 

Since  this  Is  obviously  true  for  k • N,  it  is  by  Induction  true  for  all  k ^ N. 

(Note  that,  if  N - I,  I.emme  1.1  is  identical  to  the  lemma  stated  in  Jll,  p.  17q|  , 

wlilch  is  a consequence  of  the  standard  results  for  rmrular  infinite  systems  of  alunbraic 

»iuat  ions  by  Kantorovich  ami  Krylov  in  |7,  p.  271.) 

l.et  ylx)  he  the  exact  solution  of  (1.1),  and  define  tlie  discretixation  error  function 

by  t(x)  • y(x)  - P(x)  , wliere  P(x)  is  the  linear  spline  function  approximation  to  y(x) 


-1- 


obtalnad  trom  our  numerical  method.  Denote  e(Xj^)  by  Me  state  the  following  theorem: 

Theorem  3.1.  If  y"(x)  is  Lipschitx  continuous  on  (0,1),  then  the  discretisation  error  of 
the  linear  spline  method  satisfies 


(3.1) 


- 0(h‘)  , 


- 1,  2," 


provided  the  error  of  the  starting  value,  is  of  order  h 

Proof  I By  a standard  theorem  on  Lagrange  interpolation 


(3.2) 

where 


€(x)  - - - x)e^  + (X  - x^)e^^jl  ♦ ^(x)  , 


•fiix)  - j y"(n^(x))  (x  - Xj^)  (x  - Xj^^j^) , Xj^  ^ nj^(x)  £ x^  ^ x £ x^^^ 

Since  both  y(x)  and  P(x)  satisfy  (1.1)  at  each  Icnot  x - Xj^, 

■6(8)ds  “0,  li  - 1,  2,*** 


(3.3) 


)c-l  i+1 

I / ' 

l-O  X 


/ 2 2 
i /x^-s 


This  can  be  rewritten  as 


(3.4) 


i «k  i ^i  " «it'  k - 1.  2.  •• 

i»0  ^ * 


where 

(3.5) 


k-1  l-fl 


I ! 

i-O  X, 


/ 2 2 
■'v® 


t?(s)ds  , 


and  the  w 's  are  defined  in  (2.4). 

X#  1 

Multiply  (3.4)  by  k,  difference  the  resulting  equation  for  k and  k t 1,  and  then 

divide  by  (k  + yield  the  required  error  equation 

k^i I KtI 

k 


(3.6) 

Vi- 

Jjj  Vi,i  ^i 

k - 1,  2,  •••  , 

■ 

where 

, 

■■ 

(3.7) 

k w^_^  - (k^Dw^^^^^ 

i “ 0,  !,••• ,k  , 

; 

*k+l,i  “ 

(k+l)Wk^l  k+i 

-4- 

'Vi  - \ 


k (k»l)w, 


k - 1.  2,' 


k+l,k+l 


(Note  that  (3.6)  is  not  the  only  nossible  error  eiiuation  that  can  be  deiiveil.  See,  to 
instance,  Atkinson  111  and  Benson  |21  . However  the  procedure  used  here  simplifies  the  asym('- 
totic  error  analysis  in  Section  4.) 

Equation  (3.6)  imiUies  that 


\>il  i JJVi.il  i^il  ^ l\l' 

i-0 


k - 1,  2,' 


Since  by  assumption,  6^^  • 0(h  ),  it  is  easily  shown  from  eouation  (3.4),  by  u.siivi  (..•mni.i 

2 

3.3(a)  below,  that  • 0(h  ) for  i • !,•••, K,  with  K as  defined  in  I,enma  3.2  below. 
On  the  basis  of  this  toqether  with  l.emma  3.2(b)  and  l.emna  3.3(c)  below,  we  c.rn  applv  l emm.i 
3.1  on  (3.9)  and  obtain  (3.1).  Hence  the  proof  of  Theorem  3.1  is  completed. 

Letiwa  3.2.  If  the  are  defined  by  (3.7),  then  there  exists  an  inteqer  k _ 1, 

independent  of  h and  k,  such  that 

(e>  i - i-0,  l,*"*,k-l,  k 2.  1 


“krl.k  - 


k > K , 


‘ - JJVi.il  - 


Proof  of  (a)  ; From  (3.7),  usinq  (2,4)  and  (2.5)  » wo  have 


1 1 

*k»i,o  “ nTtiiw.^,  ^ ^“'rT2 ■ V”!  ^ .2'  * 

since  the  inteqrand  is  nonneqative.  Similarly,  we  can  prove  th.it  a,  , '0  for 

k*l , i — 

i - l,”*,k-l.  By  straightforward  estimation,  it  is  easy  to  sliow  that 


" Vi.kti  ^ 


-5- 


Sine*  th«  quantity  Inside  the  square  brackets  of  the  last  expression  tends  to  (3  - 2<^)  * 0 
as  k irrereases,  there  exists  an  inteqer  K,  indei>endent  of  h and  k.  such  that 


' 0 for  k ' K. 


\»l.k  - 


Proof  of  (b) ! Since  by  usimj  (2.4)  wo  obtain 


k k 

I “k  i ■ / -7==^  ds  - j 

i-0  ' 0 .4^ 

it  can  easily  be  shown  tltat . for  k ■ 1,  2,***, 


for  k - I, 


*■  ^ 

r . ^ 1 

*k+l.i  “ ‘ ' 2 tk*l)w^ 


k+l.ktl 


1 ^ ^ 

- O riT  ■ 


By  means  of  part  (a)  and  (3.11),  the  result  of  part  (b)  immediately  follows. 

Lesia  3. 3.  If  y"(x)  is  Lipschitt  continuous  on  10,11,  tl>en  there  exist  constants 
^1'  ^'2*  ^3  ' ^ *1,  such  that 

(a)  |Pj^|  < Cj  h‘  , 


(c)  IbJ  ’ 


for  k • 1,  2,"**,  where  and  bj^  are  defined  in  (3.5)  and  (3.8),  respectively. 

Proof  of  (a)  I l.et  M,  » T**. , ly"  (x)  | . Then  by  straiqhtforward  estimation,  it  follows 

~ ■ 2 Xt  [Opli 


from  (3.5)  that 


where  C,  " t Mi* 
1 4 2 


\\\  L - 1- 


t: 


s 

I 

1 


t'ro>:)f  ot  (b)  ; Subtraction  of  (i.5)  from  (3.5)  with  k replaced  bv  k + I,  and  by  a 
change  of  the  variable  of  inte»jration,  it  is  not  difficult  to  show  that 


(3.11) 

where 


R - H - * a<2)  ^ ^(3) 


k ' 


k - 1,  2,' 


. < 11 


k-1  ’‘it2 


, * 


i “0  X , , 
1*1 


/-  ~ lv"(n^^j(s))  - y“(nj(s-h))l  (s-x^^^) (s-x^^^) ds  , 


k*l 


\ 


(2) 


, k-1  it2 

\ I / 

* “i.! 


L'^kU*"  ’Vl"” 


» X,  ,-s 
k+1 


y"(n^(s-h)) (s-x^^j) (s-x^^2)ds 


. ( J1 


X ^ i 

*0 


y"  ( njj(  s) ) (s-Xq)  (s-Xj)  ds 


I.et  !■»»  the  Lipacliitz  const-mt  for  y".  Then  by  straiqlitfoi’ward  estimation  aivl  noting 

that  hk  1,  we  obtain  ftixn  (3.13) 


Kn  - 'StI  -^2'^’  * T ■'^2  T * i "2  T-S  T-  ‘■’ 


where  C.,  - - (4l,^  ♦ IM^)  . 


Fnxif  of  (c)  : From  (3.B),  usinu  (1.12),  (3.14)  and  (2.5),  it  can  easily  be  shown  that,  for 

k - I,  2,---, 


|b^|  < k+lj_\l 


k+1, k+1 


< c ---  , 

- j .k 


where  *^3  “ T ♦ ^jl  ' 


Corollary  3.1.  if  the  assumptions  of  Theorem  3.1  are  satisfied,  then  for  anv  fixed 

X . (0,1), 


£(x)  ” 0(h  ),  £'(x)  • 0(h) 


Ti 

t 


! i 


; 


i i 


-7- 


Proof ; If  w«  express  '^(x)  which  is  defined  in  (3.2),  in  integral  fom,  wo  obtain 

X *i+l  (x.^.-s) 

v'(x)  » / {x-a)y"(s)ds  - (x-x^)  / y"(s)ds 

X.  x^ 

By  swans  of  (3.1),  Corollary  3.1  follows  ismediately  from  (3.2)  and  the  equation  resulting 
frosi  differentiating  (3.2). 


4.  An  Asymptotic  Formula  for  the  Discretisation  Error. 

In  this  section  we  obtain  an  as^TSptotic  formula  for  the  discretization  error  of  our 
numerical  solution  which  confirms  the  conjecture  bv  Noble  (101. 

Theorem  4.1.  If  y < C^lO.ll,  then  the  discretization  error  for  the  linear  spline  method 
satisf ies 

2 


(4.1) 


« Yj  y"(Xj^)  + 0(h^/»ir)  , )i  - 1,  2, 


provided  the  error  of  the  starting  value,  is  of  order  h . 

Proof ; Since  y(x)  < C^(0,1J,  it  is  not  difficult  to  show  that,  for  x t (x^, 

(4.2)  €(x)  - ^ I(X^^^  - X)(e.  - ^y»(x^))  + - TJ  Y"(x^^^))l  t ^(x)  , 


where 


'fix)  - 2 ( (x-x^)  (x-x^^j)  + — I + p(x)  , 


with 


X.  < n(x)  < X 


-,x)  - (x-x^,3  . x ,2y-(n(x^^^))  - y«(Ux.^j)))h^x-x.),  ^ ^ 


i+l 

- X , , 

- 1+1 


By  substituting  (4.2)  into  (3.3)  we  obtain 

2 


(4.3) 

Sc'  k - 1, 

i-0 

where 

)c-l  ’‘i+l 

1 . 

(4.4) 

% - - I / 

i-0  X 

-~=r  .p(s)dS 
72  2 
vx. -s 

wc  note  that 


is  converqinq  taster  than 


derive  an  error  equation  in  the  same  way  as  in  the  proof  of  converqence  without  modifvinq 


(4. a),  we  expect  the  parts  in  the  derived  error  eauation  which  corresoond  to  6 and 


in  Lenna  3.1  would  have  unbalanced  rates  of  converqence,  with  the  former  converqinq  faster 


than  the  latter.  By  applyinq  Lenxna  3.1  to  such  an  equation,  we  will  fail  to  obtain  t!ie  re 


suits  we  expect.  To  avoid  this,  we  define  ^ 


v''(x.))  and  rewrite  (4.3)  as 


Now  multiply  (4.5)  by  )t,  difference  the  resultinq  e<iuation  for  k and  k ■*  1 , then 


divide  by  (k+l)w 


to  yield  the  reouired  error  equation 


where 


k (k+l)w, 


Equation  (4.6)  implies  that 


Since  6 


it  is  easily  shown  from  equation  (4.5)  by  usinu  I.erama 


4.2(a)  below,  that  G 


K,  with  K as  defined  in  Lemma  4.1  below 


On  the  basis  of  this  toqether  with  Lemma  4. Kb)  and  Lemma  4.2(c)  below,  we  can  apt)lv  K'rain.i 


3.1  on  (4.9)  and  conclude  that 


Thus  the  proof  of  Theorem  4.1  is  completed. 


Lemma  4.1.  If  the  *154.1  *te  defined  by  (4.7),  then  there  exist  integers  K and  K ^ 1, 


independent  of  h and  k,  such  that 


“ktl.k  - 


i - 0,  1 , • • • .k-1 , k i 1 » 

k > K , 


(b)  I - I la,.  I > i 


i-0 


Vl.i'  - 4 ^ ' 


k > K , 


for  an  appropriate  C , 0 < C < 1. 


Proof  of  (a) ! By  means  of  Lemma  3.2(a),  part  (a)  follows  inmediately  from  (4.7). 

k 

Proof  of  (b)  : Since  it  would  be  very  complex  if  we  estimate  7 1 directly,  we  introduce 

k+1,1 


i-0 


*kn,i  ' 

(4.10) 


(l-C)w 


“k+l,i  “ “k+l,i  * (k+l)w 


k,i 


0<C<1,  i-0,  l,**‘,k  . 


k+1 ,k+l 

By  using  Lemma  3.2(a)  and  noting  the  nonnegativity  of  the  Wj^  from  (2.4)  , we  obtain 
(4.11)  \n,i-°'  i - 0,  i,...,k-i,  kii  , 


*k+l,k  - 


k > K , 


where  the  K is  tie  same  K as  defined  in  Lemna  3.2. 

Using  (3.10)  and  (4.11)  it  can  easily  be  verified  that 


(4.12) 


1 - I U.  - .1  > ^ 


i-0 


Vi,i'  -TTT' 


k > K 


If  we  can  show  that  for  an  appropriate  C,  0 < C < 1,  proof  is 

completed.  To  do  this,  it  is  sufficient  to  show  that,  for  i - 0,  l,***,k. 


,/k^2 


K.1,1  - '‘-''“..,1  - '"’K.i  - 1'  - “ • 


-10- 


L 


• I' 


Ity  tiAliiq  {J.4)  and  lattlnq  a • l>t,  It  l»  aaaily  verified  that,  for  1 ■ !,•••, k-3. 


, ,t)  - - (-•^  - 1)  ^ 


Sinoe  for  l-l^t^l*l,  I*  !,•••, k-3, 


l»T^c  k - rg(in)i‘*  * ia/i-c»'c  k - riTTi/tTii (in)  ^ 


If  0 is  rhoiipn  to  b«  — # it  in  obvious  from  (4.13)  that  P.  , r **  0 for  i » !,•••, k-3. 

2 krl , i — 

With  this  valuo  of  C,  It  la  also  not  difficult  to  shi>w  that  as  k incroaHos.  P.  , , 

k*  1 1 i 

(I  - 0.  k - a,  k - 1,  k)  tends  to  i-,  -( H - '2*^3  ♦ T>'2)  -ir,  4(  «.'T  - f>/i  ♦ J)  ~ , and 
4 --  1 

•jl»2  - n — , rosiwi't  t voly , and  am  therefore  qreater  than  r.eio  for  suftioir'nt  ly  lame  k. 
' >Vk 

'I'lnis,  for  i • 0,  l,-*',k,  there  exists  an  Inteaei  K ' 1,  i n<te|s>ndent  of  h and  k, 


Hiu'h  t hat 


'Vm.I  - “ktl.l' 


My  means  of  pai  t (a),  (4.11)  aiu)  (4.12),  it  follows  that 


I - !'  i“k*,,(i  1 > - ? 

i-O  ^ ’ i- 


fot  ('  • -,  and  k K,  with  K « inax(K,  R)  . 


I.emna  4.2.  If  y < O |0,1|,  then  there  exist  oonstants  O^,  i',,  i' ^ ' 0,  i luh'is'mient 


h amt  k,  sin-h  that 


(a)  lli^l  ^ • 

(o>  llr^l  :.r,  , 


for  k “ 1,  where  Rj^  and  bj^  are  defined  in  (4.4)  and  (4.8),  respectively. 

Proof  of  (a) ; By  rei^eated  integration  by  parts,  it  is  not  difficult  to  show  that,  for 
i • 0,  1 , • • • ,)t-l , 


1 1 l 

2 ^ ■ 24  / 


Using  (4.14)  we  can  rewrite  (4.4)  as 


^*i+l  x^+2s^  2 2 

f y"(x^)  (s-x^)  (s-x^^^)  ds  . 

“i  (x^-s^) 


)t  - 1,  2,- 


)c-l  ’‘i+1  x^+2s^ 

a;  > - - 24  I J — 5/2  Y"<*iH3-x^)  (8-x^^^)  ds  , 

i»0  X,  ,2  2. 

i (x^-s  ) 

,,,  )c-l  *i+l  , 

^ ‘\loL  ' 

1“U  X.  /x,  -s 
^ It 

with  p(s)  as  defined  in  (4.2). 

Let  M,  - I 1 |y"' • Then  by  straightforward  estimation, 

2 XC  10#XJ  <3  X€  lOfXJ 

and  noting  that  h)c  < 1,  we  obtain  from  (4. IS) 


l\t  1 ”2  ^ * ”3  - ^1  Tk'  k - 1,  2,-*- 


where  “ *^2  * "''s' 


Proof  of  (b)  ! Subtraction  of  (4.15)  from  (4.15),  with  k replaced  by  )t  + 1,  then  by 
straightforward  estimation,  and  noting  that  h)c  ^ 1,  it  is  not  difficult  to  show  that  for 


k - 1,  2,- 


l\ti  - 1 Ha'S  - 


< ,2„^.M3)^>5«3  5j<52j[J  . 


where  C2  ” 7M2  ♦ 6Mj. 


s 

\ I 


Proof  of  (c) ! From  (4.8),  using  (4.16),  (4.17)  ^uld  (2.5)  we  obtain,  for  )i  - 1,  2,- 

1^1  - Tirn)’!P,  ^ - \\  ^ \\U 

'M  ■ 

where  Cj  - j(C^  + 2C2) . 

5.  A Numerical  Example. 

The  linear  spline  method  was  applied  to  the  following  Abel  integral  equation: 


/ 


y^y(s)ds  - § Jo(x) 
fX  -8 


T)ie  exact  solution  is  y(x)  • cos(x)  . 

In  Table  5.1(a)  and  Table  5.1(b)  we  list  the  error  6(x)  and  €*(x)  at  )cnots  and  at 

mid-points  between  the  Knots,  respectively  on  (0,3)  for  different  stepsizes  h.  The  error 

2 

G(x)  and  G'(x)  satisfy  the  predicted  h and  h dependence,  respectively.  Note  aim'  in 

2 

Table  5.1(b)  that  the  error  e*(x)  at  the  mid-fioints  are  actually  0(h  ) although  we  have 
not  proved  that  this  will  be  the  case. 

In  Table  5.2  we  list  the  actual  error  (column  2)  for  the  linear  spline  mr'thod  and  the 
theoretical  error  estimate  (coluain  1)  cosg-Aited  from  e<)uation  (4.1)  at  Knots  on  (0,11  tor 
h - 0.01. 
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X 

> (X) 

» • U) 

h 

h ' .1 

h/'3 

h 

h/3 

h/9 

0.0 

.OOOOK  0 

.OOOOK  0 

.OOOOK  0 

. 3924E-1 

. 1 lO'lE- 1 

.43t>3E-2 

- 

0.5 

-.712».K-1 

-.7‘tt>1E-4 

-.(Ul(H.K-5 

.4345K-1 

. 145HE-1 

.4H70E-2 

1 

1.0 

-.4t.21K-  1 

-.50f.7K-4 

-.5570E-5 

.2025K-1 

.H923l':-2 

t 

1.5 

-.9744K-4 

-.'32  39E-0 

.2W<3K-2 

. 1084K-2 

. 3H23E-3 

2.0 

.2'127E-J 

. J503E-4 

-.215HE-1 

-.702 IE- 2 

• . 2 J22K“2 

2.5 

.».U'>K-3 

.704W-4 

.7724K-5 

- .4054E-1 

-.1341E-1 

-.4457K-2 

.1.0 

.7H22E-  1 

.HH73E-4 

.96fi2K-5 

-.4')5«E-1 

-.1051E-1 

-.5501E-2 

Krror  at  Mid-tnuntt)  ( h - 0 ■ 1) 


X 

e(x) 

e’(x) 

h 

h/i 

h/9 

h 

h/i 

ll/9 

0.45 

.401.2E-3 

.4 IOOE-4 

.4790E-5 

.5117E-4 

.1452E-4 

.1921E-5 

0.95 

.2341E-3 

.2f<49K-4 

,2'>07K-S 

,27<3.U>,< 

. 3352E-4 

. (HSSE-5 

1.45 

. 1410E-4 

.3312E-5 

.44  I3K-<> 

. 3075E-3 

.42»7E-4 

.4717E-5 

1 .95 

-.2070E-3 

-.20513E-4 

-.2272K-5 

. 3(,02E-  1 

.4150E-4 

.4421E-5 

-.  1775E-3 

-. 393HE-4 

-.4472K-5 

.2744E- < 

. 299  3E-4 

. UVV7K-S 

2.95  1 

-.4544E-3 

-.4rt5lK-4 

-.55H0E-5 

. 1152E-  1 

.1104E-4 

, 10'^'>E-S 
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Table  5.2 

Actual  Error  and  Theoretical  Error  Estimate 
at  Knots  with  h » 0.01 


ACTUAL 

THEORETICAL 

X 

ERROR 

ERROR  ESTIMATE 

0.1 

-.8001E-5 

-.8291E-5 

0.2 

-.7957E-5 

-.8167E-5 

0.3 

-.7790E-5 

-.7961E-5 

0.4 

-.7531E-5 

-.7676E-5 

0.5 

-.7191E-5 

-.7313E-5 

0.6 

-.6778E-5 

-.6878E-5 

0.7 

-.6296E-5 

-.6373E-5 

0.8 

-.5752E-5 

- . 5806E-5 

0.9 

-.5151E-5 

-.5180E-5 

1.0 

-.4501E-5 

-.4503E-5 
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